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The adsorption of a near critical fluid conflned in a slit pore is investigated by means of density 
functional theory and by Monte Carlo simulation for a Lennard- Jones fluid. Our work was stimulated 
by recent experiments for SFe adsorbed in a mesoporous glass which showed the striking phenomenon 
of critical depletion, i.e. the adsorption excess F flrst increases but then decreases very rapidly to 
negative values as the bulk critical temperature Tc is approached from above along near-critical 
isochores. By contrast our density functional and simulation results, for a range of strongly attractive 
wall-fluid potentials, show F monotonically increasing and eventually saturating as the temperature 
is lowered towards Tc along both the critical (p = pc) and sub-critical isochores (p < pc). Such 
behaviour results from the increasingly slow decay of the density profile away from the walls, into 
the middle of the slit, as T ^ . For p < pcwe find that in the fluid the effective bulk field, which 
is negative and which favours desorption, is insufficient to dominate the effects of the surface fields 
which favour adsorption. We compare this situation with earlier results for the lattice gas model with 
a constant (negative) bulk field where critical depletion was found. Qualitatively different behaviour 
of the density profiles and adsorption is found in simulations for intermediate and weakly attractive 
wall-fluid potentials but in no case do we observe the critical depletion found in experiments. We 
conclude that the latter cannot be accounted for by a single pore model. 

PACS numbers: 64.60.Fr, 05.70.Jk, 68.35.Rh, 68. 15. -He 

I. INTRODUCTION 

The term 'critical depletion' was introduced in connection with experiments designed to study the phenomenon 
of critical adsorption for a pure fluid at a solid substrate . When a fluid is brought to its bulk critical point in 
the presence of an attracting external wall or substrate, for example along a critical isochore, the amount adsorbed 
(adsorption) diverges as r = (T — T^jT^ 0. Tc is the critical temperature. Theory ||] attributes these divergences 
to the fact that the wall causes a perturbation of the order parameter (OP) profile miz) = p{z) — pc, where p{z) is 
the density profile, to extend over a distance ~ ^, the bulk correlation length, from the surface. Close to criticality, 
where ^\ t \~'^ {i^ is the critical exponent), the OP profile differs from its bulk value (fixed by the properties of the 
reservoir far from the substrate) over macroscopic distances and the adsorption can be a diverging quantity. Fisher 
and de Gennes |^ postulated that near criticality the OP profile should be described in terms of a universal scaling 
function, i.e. sufficiently close to Tc and for sufficiently strongly attracting walls 



m(z) = r'3A/-(0 (1.1) 



where A/" is a universal scaling function and argued that as t — > along the critical isochore (p = pc) the adsorption 
r takes on the asymptotic, universal form 

m{z)dz - r'^e - r'^'" (1-2) 

where z is the distance measured normal to the substrate, located at z = 0. (3 is the critical exponent describing the 
vanishing of the bulk OP and (3 — v ^ —0.305 for Ising magnets or fluids. In Ref. [0 measurements were performed 
for SFg adsorbed on a finely divided (colloidal) graphite adsorbent (Vulcan 2>Q). The volumetric method employed in 
these experiments measured P directly on approaching the critical point along the critical or near-critical isochores. 
Although P increased as t was reduced, consistent with ( |1.2| ), close to Tc (r ~ 5 x 10~^) the adsorption reached a 
maximum and decreased sharply taking on negative values very close to T,.. Microgravity experiments by the same 
group confirmed these results [|l],||. The critical depletion of P was attributed to the confining efi^ects of colloidal 
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particles on the near critical fluid and new experiments, designed to test this conjecture, were performed on the 
sorption of SFg in a mesoporous glass CPG-350, which comprises a rigid interconnected system of mesopores with 
a nominal pore diameter of 31nm. For the rigid porous material the measured adsorption showed a very similar 
temperature dependence to that found for the colloidal graphite adsorbent. 

These experimental results have stimulated several simulation and theoretical studies |^,^ aimed at under- 

standing the origin of the striking behaviour of F and, in particular, to answer the fundamental question as to whether 
critical depletion is a single pore phenomenon, as is for example capillary-condensation which arises from the com- 
bined effects of wall-fluid forces and finite-size and which manifests itself in simple confining geometries. The grand 
canonical Monte Carlo simulations |^,|| of a Lennard-Jones fluid confined between two structureless, attractive planar 
walls did indicate that the average, local density in the middle of this slit-pore could fall below pc under near-critical 



conditions thereby leading to depletion. However, it has been shown recently that the depletion found in [4|j^] is an 
artifact of the simulation procedure |27). The theoretical studies [|| of the simplest model of a confined fluid, namely 
the lattice gas subject to identical surface flelds located at the two walls, revealed a physical mechanism which could 
cause a dramatic decrease of the adsorption F on approaching Tc along near-critical isochores. Consider the situation 
when the density of the bulk reservoir is slightly lower than critical. Transcription of the lattice gas into the Ising 
model sets the chemical potential difference 

^^l = ^l{p,T)~^l{p,,T)^2H, (1.3) 

where p{pc,T) is the chemical potential on the critical isochore and H is the bulk magnetic field. Thus for p < pc, 
A/i ~ 2H < and the bulk field favours the dilute (gas) phase. If the surface fields are sufficiently attractive that 
they favour adsorption of the dense (liquid) phase then one has a competition between bulk and surface fields which 
influences the shape of the OP profile and hence the behaviour of the adsorption in the slit-pore. When the pore 
is large and when the bulk correlation length is much smaller than the width of the pore L, then the fluid in the 
middle part of the pore should behave as a bulk fluid. For weak H , the bulk OP (magnetization) in the critical region 
behaves as 

rub^Hx^ Ht-"' (1.4) 

where x is the susceptibility. Near the walls, on the other hand, the fluid should behave as in the semi-infinite near 
critical system subject to the surface field. 'Bulk' and 'surface' fields give diverging contributions, but of opposite 
sign, to the adsorption F for t 0. For large pore widths and for t ^ tq, where tq is defined by ^(tq) ~ L, the 
adsorption can be approximated, for < 0, by 

r ^ AiT^-" - A2\ H It-'L, (1.5) 

where Ai and A2 are positive amplitudes. Since 7 > (ly — P) the second term always dominates for T sufficiently close 
to Tc provided | | is sufficiently strong and, as a consequence, depletion of F will occur. For even smaller r, such 
that r ^ To, the adsorption saturates at a value which depends on H and on the surface field. The lattice gas model 
of a single pore predicts critical depletion for p < Pc only. For the case of the critical isochore p = pc, the bulk field 
H = and as t is reduced F first increases monotonically following the Fisher-de Gennes power law and eventually 
saturates taking on a positive value at r = (T = Tc) 

This scenario was confirmed by explicit mean-field lattice gas calculations |^ and for two dimensional Ising films, 
with bulk and surface fields of opposite sign, by density matrix renormalization group calculations j^]. In the latter 
case it was shown that the near-critical fluctuations can lead to an even richer variation of F(r). For certain values 
oi H , in addition to the maximum, a minimum of F appears where the correlation length approaches the pore width 
and the competition between the effect of adsorbing walls and the large susceptibility of the central region (favouring 
the dilute phase) becomes particularly strong. For weak H < the results for F(r) obtained in the lattice gas (Ising) 
model of a single pore look very similar to those measured in the experiments of Ref. ||] and (2| . In the experiments 
for SFe adsorbed in the controlled pore glass the actual isochores correspond to densities lower than critical, i.e. 
p/pc = 0.995 and 0.999, so that < 0. 

It is tempting then to argue that since A/i < there is an effective bulk field H < which competes with the 
surface fields to give rise to critical depletion. Although this is a plausible explanation of the observed phenomenon 
it does not take into account the actual situation in a fluid. For example, if the bulk density p is fixed according 



to the experimental condition of the fluid reservoir, A/i as defined by (1.3) will vary as T approaches Tc- The 
corresponding bulk magnetic field will vary in the same way. In the present paper we show that taking into account 
the temperature dependence of Ap has a profound effect on the behaviour of confined fluids near Tc- In particular we 
flnd that under the experimental conditions of Ref. ||l|,|[ a simple fluid confined in a single slit pore should not exhibit 
critical depletion. Rather the adsorption should increase monotonically as r — > 0. This implies that an explanation 
of the experimental observations is still lacking. 
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Our paper is organized as follows. In Sec.2 we reconsider the physical mechanism which leads to critical depletion of 
adsorption in the case of the lattice gas model of a pore considered in Ref. |^ and give a heuristic argument as to why 
this phenomenon should not be expected for real fluids when the reservoir density is fixed according to experimental 
conditions. Our argument is supported by explicit calculations of T using the density functional approach and by 



grand canonical Monte Carlo simulations of the Lennard-Jones fluid in a slit-pore. In Sec. Ill we report density 
functional results obtained from a square gradient approximation to the free energy functional and short-ranged 
(contact) wall-fluid potentials. Both classical and non-classical bulk free-energy densities are employed and and in the 
classical case we investigate two forms of the free energy density, namely the Landau model free energy and the free 
energy of the Lennard-Jones fluid as obtained in Ref. |^ from an accurate fit to simulation data. Sec. 4 describes the 
computer simulations of the density profiles and adsorption of the Lennard-Jones fluid on the critical isochore and 
for a sub-critical isochore p < Pc- Results are presented for various strengths of the 4-10 and 3-9 wall-fluid potentials. 
We conclude in Sec. 5 with a discussion of the relevance of our findings for the experiments described earlier. 



II. HEURISTIC ARGUMENT. 

Here we reconsider the scaling argument Q that predicts critical depletion in the lattice gas model of a single pore 
and modify it to incoroprate two features that are relevant for the case of fluids. 

First, in order to mimic the experimental situation more closely one should account for the fact that Afi = 
p{p,T) — p,{pc,T) will vary as the temperature T approaches Tc at constant p. This implies that the corres pond ing 
bulk magnetic field H should also vary with T and therefore that the second term of the approximate formula ( |l.5[ ) for 
the adsorption T might have a different r dependence. This in turn may affect the result of the competition between 
the bulk and surface fields. 

Second the lattice gas model considered in Ref. |^ has an exact particle-hole symmetry, which corresponds to the 
trivial symmetry under reversal of the field H in the equivalent Ising model. For real fluids such symmetry is only 
approximate. It is well established that the reduced symmetry of fluids leads to scaling field mixing close to the 
critical point ||l0|| . To linear order in r and in p[p^ T) — /i^, where Pc = m(Pcj 7c)j the scaling fields have the form 

UH ^ p{p,T) - Pc ~ ClT (2.1) 



T + C2{p{p,T) - Pc), 



(2.2) 



where the parameters ci and C2 are system-dependent (non-universal) quantities controlling the degree of field mixing, 
ci is identified as the limiting critical slope of the coexistence curve, i.e. ci/Tc — dp'^°'^^{T)/dT. In order to 

account for the asymmetry of a real fluid near its critical point one should identify the bulk field 2H with uh rather 
than with Ap in the scaling analysis. 

The temperature dependence of the bulk field H depends on the particular equation of state. Consider first the the 
simplest possibility, i.e. the classical equation of state in the critical region given by retaining only the leading terms of 
the expanded van der Waals (vdW) equation of state. In terms of reduced temperature r and density r = {p — Pc)/Pc 
the vdW equation of state reads 



Ap* 



-6r 



I In 



1 - r/2 
1 + r 



+ 4(1-Kr) 



1 



1 - r/2 



1 



(2.3) 



where Ap* — Ap/PcVc, Pc is the critical pressure and Vc is the critical volume per molecule. The leading order 
behaviour of this equation in the near critical region is 



A^* = 6rr 



(2.4) 



where we have ignored terms O(r^) and O(rr^) and higher. Note that Eq.(2.4) exhibits particle - hole symmetry in 
that Ap* along an isotherm is antisymmetric with respect to the critical isochore 

Ap*{~r, t) = -Ap*{r, r). (2.5) 

Moreover, for this case there is no scaling field mixing and uh ~ Ap* whose magnitude decreases linearly in r as Tc 
is approached at constant r. 

In order to analyse the influence of a r-dependent bulk field on the behaviour of the adsorption T(t ) we reconsider 
the approximate formula (|l.5|). For the expanded vdW equation of state (2^) and r < {p < pc) Eq.(|l.5|) becomes 
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As 



(2.6) 



where ^3 = (3/4)^2, A4, = 3A2 and for consistency with the vdW approach the critical exponents should take on 
their classical values, (3 = v = 1/2, 7 = 1 and the first (critical adsorption) term diverges as -Inr. Apa rt fr om the 
additional temperature independent term, Eq. (|2.6| ) has the same form as for a constant bulk field (Eq. (|l.5| )). The 
additional term does not affect the shape of the curve r(T). It simply shifts r(T) as a whole towards negative values 
and for large widths of the pore L and (or) relatively large |r| it could drive T negative sufficiently far from the critical 
point (large r). Closer to Tc the temperature dependent terms dominate and the analysis of T{t) as a function of H 
and L performed in Ref. |^ for constant bulk field H, goes through with \H\ replaced by |rp. Following Ref. [|j we 
rewrite Eq.(^.6|) in the form 



r(r) 



Ax 



A4\r\L, 



(2.7) 



where 



Tr = (Aslrl^L) 



3rU/(A-iy) 



(2.8) 



and we have used the exponent relation pl[ 7 = A — /3 to introduce the gap exponent A. Once again it is understood 
that the exponents take their classical values. In the region of validity of approximation (1.5), i.e. for 1 S> t ^ rg, 
with ^(tq) ~ L three different ranges of with qualitatively different behaviours of T{t) can be distinguished 

1. Tr ^ To 



In this case the first term in square brackets in (2.7) dominates in the whole region of validity of this approxima- 
tion since t/t^ ^ 1 and v — A is negative. Hence, in this region the adsorption should increase monotonically 
as r ^ 0. 

2. > 1 

This condition is equivalent to |rpL 1. In this case t/t^ ^ 1 throughout the critical region and the second 



term in square brackets in (2.7) dominates. Hence, the adsorption is negative and desorption takes place despite 
the presence of adsorbing walls. 

3. To < Tr < 1 

For a given pore width L , the second term in square brackets in (^^) dominates so long as t < t^ and then 
r(T) is negative. As t inc reases, r(T) reaches a maximum for t w t^.. Finally for t ^ t^ the constant term 
in square brackets in (2.7) dominates over the second term, and for such temperatures the usual Fisher - de 
Gennes type of adsorption should occur. 

We now consider values of parameters appropriate to the experiments of Ref. 0,^. Assume that L/cr ~ 10^, is of 
the size of the nominal pore diameter ( a is the molecular size) of the mesoporous glass used as the adsorbent. For 
r — —0.001, —0.005, corresponding to the two near-critical isochores along which r(T) was measured, t^ ^ 10^^, 10^^. 
For this value of L, tq ~ 10""* and t^ ^ tq. Then according to the above discussion the adsorption F should increase 
monotonically as T is lowered towards Tc following these two isochores. Even for densities that deviate more strongly 
from pc e.g. r = —0.01, t^ ~ lO^"' ~ tq and the condition (3.) for the occurence of depletion of F might still not 
be satisfied. It is also very likely that for these values of r and L the adsorption is positive for t « To and hence for 
T < To it should saturate at a positive value. 

It is important to contrast this constant r (density) scenario with the constant H lattice gas described in Ref. 
There depletion was observed for fields H in the range -10""^ to -1.5 x 10~^. These were sufficiently strong to drive F 
negative for t > tq, i.e. while ^ was smaller than L. F then saturated at a negative value for t < tq. In the present 
case, even for values of r as negative as —0.01, the effective field might not be strong enough to drive F negative 
before ^ ~ L and then F would saturate at a positive value. 

O ur analysis so far has been based on (2.4). Consider now equations of state which do not incorporate the symmetry 
(2.5) in the fi — T plane. For systems described by such equations of state the 'true' scaling fields are now uh and 
Ur and by analogy with bulk the 'true' OP m(z) which satisfies the scaling relations (1.1) and (1.2) is not p — Pc 
but rather the linear combination of the number and entropy densities {p — pc) — €2(3 — Sc). The entropy term in the 
OP does not change the leading asymptotic behaviour of the adsorption for the semi-infinite system; it gives rise to a 
"correction term" to Eq.(p^) of order 7-1-"-'^^ where a is the specific heat critical exponent. In the scaling analysis 
of F for the confined system, the scaling field H should be now replaced by uh- In order to see if this can change the 
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behaviour of T we first consider equations of state that are hnear in r, as was the case in the vdW equation of state. 
For all such equations the 'mixed' scaling field uh reduces to 2H defined by ( |l.3| ). This is due to the fact that the 
chemical potential on the critical isochore and the chemical potential at coexistence have the same limiting slope at 

Tc, i.e. ci/Tc = lim^^y- (d^™^^/dT) = limj,^^+ {dfi{p,T)/dT)p^ so that 

Uh = T) - fic- ciT ^ Afi + ^{pc, T) ~ pc - ciT = A/i = 2H. (2.9) 

Thus, for the classical equations of state which are linear in temperature the reduced symmetry of the fluid does not 
influence the temperature behaviour of the bulk flcld. 2H — A/i differs from the leading order behaviour only by 
terms higher in r. For example, in the case of the vdW equation of state these are of order r"* and rr^. 
If the equation of state is not linear in temperature then 

UH = Ap + p{pc, T)-pc- ciT ^Ap + a2T^ + 0{t^), (2.10) 

where 02 is a constant coefficient. Thus, the temperature dependence of the scaled field uh, which now plays the role 
of the bulk field H, differs from that of Ap, but only by terms higher order in r. 

We conclude that differences arising from the reduced symmetry of the fluid, i.e. 'mixed' scaling flelds, are not 
important for the behaviour of the adsorptiom ^The presence of higher order terms in r and r do not change the 



conclusions of our analysis performed using (2.4). Thus our predictions of no depletion of adsorption along near- 
critical isochores should be valid for all classical equations of state. 

Of course real fluids are non-classical. Our argument can be extended using the fact that near criticality real fluids 
should obey the scaled equation of state jl^ 

Ap = r\r\^'^Doh{T/\r\^/'^) (2.11) 

where S is the critical exponent and Dq is an amplitude for the power law on the critical isotherm, and h{x) is a 
scaling function. Note that Ap = p{p,T) — p^°'^^{T) for r < and 5 = 1 -\- "(1(3- Although there is no a priori 
theoretical expression for h{x) the scaling function should satisfy several conditions following from requirements of 
thermodynamic stability and analyticity of the chemical potential. Thus, h{x) should be analytic in its range of 
definition —1 < x < 00, equal to at a; = —1, the coexistence curve, and possess an (asymptotic) series expansion 
near a; = 00 (the critical isochore) of the form 

00 

/i(.t) = ^ry„x^(^+i-2„)^ (2.12) 

n=l 

For small values of x, h{x) should have an expansion of the form 

h{x) = l + ^Kx''. (2.13) 



The leading temperature dependence of Ap on the near-critical isochores is given by the first term in expansion ( 2.12| ) 
and the leading r dependence is given by the first term in expansion ( 2.13| ), i.e. 



A/i - ?7irr'^ + L>or|r|*"\ (2.14) 



where Dq and 771 are amplitudes. For classical exponents this expression is consisten t w ith Eq.(2^). Using this 



form for the bulk field 2H ~ Ap we can repeat the analysis performed above. Eq.(2.7) remains valid but now 
Tr - (|r|*L)i/('^-''). For real fluids 5 - 4.78 which means that the values of are even smaller for a given r than in 
the classical case. This implies that the effective bulk field is very weak for the conditions of the experiment and thus 
depletion of the adsorption F should not occur. Rather saturation of F at positive values should be expected. 



III. RESULTS FROM DENSITY FUNCTIONAL THEORY 



In this section we report the results of density functional calculations for the adsorption F of a near-critical, simple 
fluid confined in a slit-like pore. These results provide an explicit test of the heuristic ideas given above. 

Specifically, we consider a fiuid confined between two identical parallel adsorbing walls located at z = and z = L 
and infinite in the x and y directions. The system is in contact with a bulk reservoir at fixed temperature T and 
chemical potential p. The equilibrium profile is obtained by minimizing the grand potential functional [T^ 
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n[p]=T[p\- j dr{fi-V{r))p{r), 



where ^(r) is the total wall-fluid external potential 

V{r) = V{z) = U{z) + U{L - z) 



(3.1) 



(3.2) 



and U{z) is the solid-fluid potential due to a single wall. The equilibrium density profile p(r) = p{z) corresponds 
to the minimum of ^[p]. We choose the simplest form for f2 based on the square gradient approximation to the 
intrinsic free energy functional J-[p\ and model the wall-fluid contribution by a term $s which depends only on the 
fluid density at contact i.e. on p{Q) = p{L). In this approximation the grand-potential excess per unit area is the 
following functional Q : 



l[p\ = 




D (dp 

dz 



$3 



(3.3) 



Here tpi^p) = w(p) -t- P is the excess grand-potential density, i.e. Lo{p) = f{p) — up is the grand potential density, P is 
the pressure and f{p) is the Helmholtz free energy density of a homogenous fluid of density p. For T < ipip) has 
two minima corresponding to the two distinct bulk phases. At bulk coexistence both minima are equal to zero. The 
coefficient D is related to the second moment of the direct correlation function [|l3| but for simplicity we choose it to 
be density independent. The wall-fluid term has the form 



Up^iO) + p\L)) - eUpiO) + P{L)). 



(3.4) 



The first term, with c > 0, represents a reduction of attractive pair interactions between fluid particles at the surface 
arising from exclusion of the fluid by a wall. The second term with > measures the strength of the attractive 
wall potential. Symmetry of the wall-fluid potential dictates that p(0) = p{L) and dp/ dz = at z = L/2. 

It is well known that functionals of this type cannot incorporate short-ranged correlations and hence cannot acount 
for oscillations of the density proflle which occur near the walls |^. However, they should capture the main features 
of critical adsorption in large pores, as this phenomenon is dominated by the behaviour of the proflle far from the 
walls. Indeed they were succesfuUy employed by Marini Bettolo Marconi ||l^ in a pioneering study of the effects of 
flnite size on critica l ad sorption. 

Minimization of (|3.3|) yields the following equation for the density profile p{z): 



D 



d^p{z) dip 



dz^ 



with boundary condition at the wall z = 0, 



D 



dp{z) 
dz 



dp{z) 



cp(0) - £^ 



z=0 



Eq. (p.SD has a first integral 



D 



dp{z) 
dz 



= ^ip) + F, 



(3.5) 



(3.6) 



(3.7) 



where P is a constant of integration, independent of whose value depends on T,L and p. The function F{L), which 
vanishes as L ^ oo, can be identified with the solvation force between the walls |p^ , |l6| i.e. F{L) = —2{dj/dL)T,^, 
where 7 is the equilibrium value of 7[p]. It may be determined from the equations 



L = {2Dy/^sgnipiO) - p{L/2)) 



p(o) 



dp 



p(L/2) mp)~ML/2))Y/^ 



and 



F 



'^Ap[L/2), 



(3.8) 



(3.9) 



both of which follow from Eq. (|3.7|) along with 
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^'(^(O)) - Hp{L/2)) = ^(cp(0) - e^f (3.10) 



which foUows from the boundary condition ( ) . The key quantity of this study, the Gibbs adsorption Tg (coverage) 
is defined as 

Tg - / {p{z)~pb)dz (3.11) 
Jo 

with pb the density of the bulk fluid at chemical potential p and temperature T . Tg satisfies the Gibbs adsorption 
equation 

Vg - -2[d-i{L)/dp)T. (3.12) 
which gives, using Eq.(3.3), the following expression for Tg 

= {2Df'\.gn{p{0) - p{L/2)) '^.'^'u^n/oin/^ - (^.13) 

Jp{L/2) mp) -'^(p{L/2))]y^ 



T 



G 



In order to test our predictions from Sec. 2 we chose three different models for tpip) ^^^d calculated the adsorption 
as a function of temperature on approaching Tc from above along near-critical isochores, i.e. for fixed p < pc- 

(a) Landau model free energy. 

In this case we expand the grand potential density uj{p) and the pressure P — —uj{ph) about the critical density pc- 
In terms of reduced variables r — {p — pc) / pc and j't, = (pb — pc)/ pc the dimensionless excess grand potential is 

r{r) = y (r^ - rl) + ^ (r^ - r^) - (r - r,)^p\ (3.14) 

where tp* = ^{p)/kBTcPc, {kBTd Pc)a* = a = {dp/dp)T at p = pc, [kBTd pl)h* =b= \{d'^p/dp^)T a.t p = Pc and 

^p* = {p-p{pc.T))/kBTc ||. 

Such a choice for 'ip{p) corresponds to the simplest mean- field or Landau description of a model fluid exhibiting 
particle-hole symmetry. For the special case A/i = (on the critical isochore pb — Pc) the integrals (3.8) and ( 3.13| ) 



can be performed explici tly in terms of the incomplete Jacobi elliptic integral of the first kind [M. For A/i < 0, i.e. 



for Pb < Pc, the relation (3^) between the order parameter at the midpoint r„i and the wall separation L can also be 
expressed in terms of the elliptic integral of the first kind (see Appendix). 

For a given value of the separation L between the walls and for fixed bulk density pb, we determined the reduced 
densities and Tu, and F{L) a t va rious temperatures, corresponding to r between 0.1 and 0, using a graphical 
construction [|l6| along with Eq.(|3.8|). For a given L there is only one solution in this range of temperatures. At 



each temperature we calculate numerically the integral (3.13) for Tg using the Romb erg metho d. A s a check of the 



accuracy we calculated Tg for r;, = using analytical expressions for the integrals (3.8) and (3.13)in terms of the 
incomplete Jacobi elliptic integral of the first kind. 

We performed our calculations using parameters a* and b* obtained by fitting the critical temperature and the 
critical density of SFg, i.e. pc = 3.05 nm~^, Tc = 318. 7K using a generalized vdW equation of state (see Ref. ||l^ for 
details). We find a* = 2.764r and b* = 0.113. The value of D was taken from Ref. (l|], giving D* = 1.094. The waU 
fields c and e^, were treated as independent parameters which were varied in order to examine the influence of the 
strength of the wall potential on the behaviour of Tg (t) . 

We studied systems with different wall separations L ranging between 25 and 200 nm. For each L and fixed wall 
fields c and we calculated Tg{t) along the critical isochore and for several near-critical isochores, i.e. for fixed rb 
between and —0.1. 

The results for L = 100 nm, c* = c/D — 0.5 nm~^ and ejj, = e^/ipcD) = 1 nm^^ and several values of rb are 
shown in Fig.lA. Similar behaviour of the adsorption is found for the other values of L we considered. Note that the 
effective hard-sphere diameter resulting from fitting the SFg data is 0.433 nm [ p^ . In order to allow for a more direct 
comparison with the analysis of Sec. 2 and with the results from the lattice gas (Ising) model of Ref. Q we plot in 
Fig. 1 the quantity 



L 

dz{p{z) ~ Pc) = Tg + PcUL (3.15) 
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rather than Gibbs adsorption Tg- Our choice for the coefficients c* and corresponds to walls that attract fluid 
rather strongly. On the critical isochore, i.e for rb = 0, the reduced contact density r„ ~ 0.57 at T = Tc and does 
not change very much as the temperature is increased. Our results show no depletion; for all choices of the the bulk 
density the curves ^dT) increase monotonically with decreasing r. The adsorption first increases as T approaches 
Tc from above and then saturates sufficiently close to Tc- Except for = —0.1, the value at saturation is positive. 
Moreover for the isochores which are the nearest to the critical one, i.e. for < < —0.01, the values of Tdr) are 
very close for r < 10"'* and the adsorption saturates at r « 10~^ corresponding to ^ ~ L, which agrees with the result 
found from the Ising model We have examined the behaviour of Tdr) for several choices of the coefficients c* 
and £*, but no qualitative differences in behaviour have been found. When the strength of the wall-fluid interaction 
decreases, the degree of adsorption becomes smaller. 

(b) Lennard-Jones fluid free energy. 

As a second model we chose f{p) to be the Helmholtz free energy density of the Lennard-Jones (LJ) fluid as given 
by the empirical, modified Benedict- Webb-Rubin (MBWR) equation of state 

np*)^P*A*{p*)^P* (E^+E&.gJ +p*A*, (3.16) 

\i=l 1=1 / 

where f*{p*) = f{p)<T'^/e and p* = pa^. e is the LJ well depth and a is the LJ atomic diameter. A* As the ideal gas 
term, and bi are nonlinear functions of temperature and Gi are functions of the density p* (see Ref. ^). The MBWR 
equation is a classical equation of state obtained by fitting to simulation data. It describes the near-critical region 
of a LJ fluid quite accurately. Since it is non- linear in temperature it does not possess the particle-hole symmetry of 
&■ 

In order to calculate the adsorption in this model we proceed along the lines described above for the the Landau 
model but now we perform all the integrals numerically. We fit the parameters to the critical temperature and critical 
density of SFg and the parameter D is obtained from the relations a = {dp,/ dp)Tc sd p — pc and 

a/i5 = e„-V2-, (3.17) 

with I' = 1, the classical value. the correlation-length amplitude, is set equal to the experimental value 0.2 nm for 
SFg. We find D/{a^e) = 0.364 nm^ and a = 0.467 nm. 

Fig. IB shows the results for F^ for the wall separation L = SOcr and ri, between and —0.15. The wall fields 
are chosen to be equal to c* = c/D = 0.5 nm~* and = ew/{pcD) = 0.75 nm~*, which yields a contact density 
Tw ~ 0.51 at the critical point. The shapes of the curves are very similar to those for the Landau model. The absolute 
values of Fc are smaller as L is much smaller than the value used for the Landau model. Note that for this wall 
separation the density of the bulk reservoir must be further removed from pc in order to shift the whole adsorption 
curve below zero; for = —0.15, Fc still saturates at a positive value. 
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FIG. 1. Adsorption Tc{t) (in units of nm~'^) defined by Eq. ( p.Iq ) as a function of the deviation from the bull; critical 
temperature, calculated in three different models (A) Landau model for wall separation L — 100 nm, surface field = 1 
nm~^ and surface coupling c* = 0.5 nm~^. The curves correspond to different values of the reduced bulk reservoir density 
Vb = {p — pc)/pc- From the top to the bottom: rj, = 0, —0.005, —0.01, —0.03, —0.05, —0.1; (B) model with the free energy 
density of the LJ fluid for wall separation L = 50a, where a is the LJ atomic diameter, surface field ej^, = 0.75 nm~^ and 
surface coupling c* = 0.5 nm^^. The curves correspond to the following values of rt: rt = 0, —0.005, —0.01, —0.1, —0.15; (C) 
Fisk-Widom model for the same L, e^, c* and rt as in the Landau model. 



(c) Fisk-Widom free energy. 

In order to incorporate non-classical critical exponents we employ the simp lest possible approach that goes b eyond 



mean field, i.e. the Fisk - Widom functional ||17|. This has the form of (3.3) with the grand potential density (3.14) 
replaced by 

r{r) = y (r^ - rg) + ^{r'+' - 4+') - [r - n)Af,*, (3.18) 

where ^j* = ^{p)/Pc, [Pc/ pl)a* =a = {dfi/dph at p = p„ [PJ pi+^)b* = h and A/^* ^ {p - p{p,,T))p,/ P,. The 
parameter a now vanishes as the 'exact' inverse compressibility: a ~ r'' and the dimensionless quantities are defined 
in a slightly different way from those in the Landau model - see below ( |3.14 ). Following Ref. we invoke rational 



approximants for the critical exponents, i.e. 7 = 4/3, p = 2/3, /3 = 1/3 and (5 = 5. In this approximation the specific 
heat exponent a and the correlation function exponent rj are equal to zero and the coefficient D can then be treated 
as a constant, whereas in reality D diverges as t~'^'^ . We define D* = Dp^/Pc, which has dimension length^. 

As for the Landau model, the integrals in the formulae for L and Tq can be expressed in terms of elliptic integrals 
of the first kind for the special case of the critical isochore = (see Appendix). We use these formulae to calculate 
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Fc along the critical isochore. For A/i ^ we evaluate the integral for F numerically. The parameters D*, a* and 
b* were obtain ed, following Ref. ||l^, from the Missoni-Levelt Sengers-Green nonclassical equation of state and from 
relation ( 3.17 ). These were fitted, as in the previous models, to the experimental values of Tc, pc and for SFg. The 
critical pressure Pc — 3.7605 MPa. 

Fig.lC shows our results for Fc calculated for the same wall separation L = 100 nm, the same values of the wall 
fields c* — 0.5 nm~^ and = 1 nm^^ and the same choices of the reduced bulk densities r?, as the results obtained 
from the Landau model and shown in Fig.lA. Again, there is no depletion and the overall form of the results for 
the adsorption is similar to that for the Landau model and for the model with the LJ fluid free energy density. The 
adsorption is stronger since j3 — v = —(1/3) for this model rather than (logarithmic increase) for mean field. After 
saturation has set in, i.e. for r < 2 x 10^^, Fc is almost independent of r^. From the plot of the midpoint density 
as a function of temperature (Fig. 2) we see that the temperature dependence of the adsorption along the isochores 
mimics the temperature dependence of r^. takes on a value equal to Th for large r and then increases as t — > 0. 
This behaviour is completely different from the case of the of the lattice-gas (Ising) model with constant negative 
bulk field in which the magnetisation becomes more negative as Tc is approached from above. 
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FIG. 2. The reduced midpoint density r™ = {p{L/2) — pc) /pc as a function of the deviation from the bulk critical temperature 
calculated in the Fisk-Widom model for the same parameters as in Fig. IC. 

IV. SIMULATION STUDIES OF A LENNARD-JONES FLUID 
A. Computational details 

We have performed Monte Carlo simulations of a simple one component fluid, interacting via an interparticle 
potential of the Lennard-Jones (LJ) form: 



ULjir) 



4e 



(7)"-(7)' 



(4.1) 



Here, as in Sec. HI, e measures the well depth of the potential, while a sets the length scale. As is customary in 



simulations of systems whose interparticle potential decays rapidly with particle separation, the LJ potential was 
truncated in order to reduce the computational effort. In accordance with many previous studies of the LJ system, 
the cutoff radius was chosen to be rc — 2.5cr, and the potential was left unshifted. No corrections were applied to 
account for the effects of the truncation. 

The simulations were performed within the grand canonical ensemble [ p^ , pO| , permitting fluctuations in the total 
particle number N. Two distinct geometries were studied: 

(A) A fully periodic cubic system of volume V — L'^. 



(B) A slit-pore geometry in which the fluid is confined to a cuboidal simulation cell of dimensions x Ly x L (with 
Lx = Ly) having structureless hard walls in the planes z = and z = L, and periodic boundary conditions at 
the cell boundaries in the x and y directions parallel to the walls. 
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Consider the behaviour of the configuration averaged local number density p(r) in these systems. In geometry (A), 
translational invariance ensures that p{r) is independent of the position vector r{x, y, z), and one has simply p(r) = p, 
i.e. the configuration averaged number density. In the absence of finite-size effects, p is completely determined by the 
imposed values of the chemical potential p and the temperature T. By contrast, in geometry (B) the presence of the 
walls at z = and z = L break the translational symmetry in the z direction giving rise to a onc-dimensional density 
profile p{r) = p{z) representing the configuration averaged local number density at a given z. The precise form of this 
profile depends not only on p and T, but also on the details of the fluid-wall interaction. In this work, we assume 
that fluid particles interact with a single wall via a long ranged potential having one of either two forms: 



U4-io{z) = 4e/4-i 



(4.2a) 



C/3_9(z) =4e/3- 



2 /CT 

15 v7 



(4.2b) 



where / is a parameter that tunes the strength of the fluid-wal l int eractions relative to those of the fluid interparticle 
interactions. The total wall-fluid potential is then given by (3.2). We note that U4-io{z) models a wall which is 
assumed to comprise a single plane of LJ particles, while U^-g models a wall that fills the half space |^l|. Since both 
of these wall potentials decay considerably less rapidly with increasing separation than the LJ interparticle potential 
of eq 



4.1 



no potential truncation was applied. 



It is instructive to compare the forms of the two types of fluid-wall potentials (eq. 4.2) both with one another, 
and with the Lennard- Jones (6-12) interparticle potential (eq. 44). This comparison is made in fig. || for the case 
4e = / = 1. The relative range of the two wall potentials is exposed by the inset which shows the result of first scaling 
the well depth of the 4-10 potential to equal that of the 3-9 potential, and then translating along the abscissa until 
both minima coincide. In this representation, one sees that the 3-9 potential is both longer ranged and exhibits a 
broader minimum than the 4-10 potential — a fact which will aid the interpretation of the simulation results presented 
below. 
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FIG. 3. Comparison of the potentials of eqs. and 4.S for the case 4e = / = 1. The inset shows the result of scaling the 
well depth of the 4-10 potential to equal that of the 3-9 potential, and translating along the abscissa until the minima coincide. 



B. Determining the critical isochore 

As far as practicable, our simulation strategy has been to try to mimic the experimental adsorption studies of rcf |^ 
in which r(r) was measured along the bulk critical isochore as Tc was approached from above. A prerequisite in this 
regard is an accurate knowledge of the locus of the bulk critical isochore. Within our grand canonical simulation 
framework, this is specified by the function p{pc, T). For the LJ fluid (with = 2.5cr), p^ is known accurately from a 
previous finite-size scaling (FSS) study, as are estimates for p^. and p^ ]. Hitherto, however, no accurate estimates 
for the critical isochore itself have been reported. Accordingly a new set of simulations were performed to determine 
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/i(/9c, T) for a range of super-critical temperatures. These simulations were carried out using simulation geometry (A) 
described above. 

The procedure adopted for estimating fi(pc,T) is detailed below and involves determining, for a given temperature 
T, that value of the chemical potential for which the measured density matches the known critical point value. 
Unfortunately this task is complicated by finite-size effects. Well away from the critical point the correlation length f 
is small, and provided the linear dimension of the periodic simulation box L ^ ^, the finite-size function fi^^pc, T) will 
provide a reliable estimate for the bulk isochore /i(pc, T). As the critical point is approached, however, the correlation 
length grows until it is comparable to, or greater than the system size L. In this regime, one expects that estimates 
of ^l{PctT) will deviate systematically from the limiting bulk form. 

In principle, finite-size scaling (FSS) methods can be employed to obtain estimates for bulk quantities from sim- 
ulations of finite size . Unfortunately their application near the critical point is rather less straightforward than 
at the critical point. The difficulties stem from crossover effects associated with non-zero values of the two relevant 
scaling fields |10|, uh and u^, which control deviations of the number density and energy density from their critical 
point values-see Eqs. (2.1) and (2.2). Small but finite values of these fields result in a large but finite correlation 
length ^ which, owing to computational restrictions on the range of accessible system sizes, renders the limit L 3> C 
(in practical terms) unattainable. One is therefore forced to attempt to extrapolate to the thermodynamic limit using 
data from system sizes for which L ^ ^. However, such an extrapolation is fraught with complications since it requires 
prior knowledge of the universal scaling functions (and associated non-universal amplitudes) controlling the crossover 
to the thermodynamic limit as both \ur\ and \uh\ are increased. To our knowledge, accurate forms for these scaling 
functions are not available. 

In view of these difiiculties, we have not attempted a full FSS analysis of the critical isochore. Instead we have 
simply determined ij,l{pc,T) for the largest accessible system size and used this as our estimate for the bulk function 
fi{pc,T). Notwithstanding the lack of a FSS analysis, there are formal grounds for believing that for the special 
case of the critical isochore, finite-size effects are smaller than on any other near critical isochore. To see this one 
must consider the effect on observables of finite values of the scaling field \uh\- Specifically, we focus on the ordering 
operator conjugate to uh, which is given by ~ (p — C2s), representing some particular linear combination of the 
number and entropy densities p5| ] . 

The dependence of SM = A4 — Aic on \uh\ differs between the thermodynamic and FSS limits. For the former 
case one has 5M ~ u'h^ , while for the latter, 5M.L ~ uhL^^^ ■ It follows that for the special case uh = 0, estimates 
of SA4l will agree with the bulk value of SA4 for all L, i.e. exhibit no finite-size dependence. For non-zero \uh\, 
however, there will be a finite-size error SA4l — SA4. Specifically, for a given L one expects that as \uh\ is increased 
from zero, SMl — SA4 initially increases like but slows with increasing \uh\, reaching a maximum at some 

\uh\ ~ L~^^/^. Thereafter, further increase in \uh\ lead to a decrease in the magnitude of 5M.L — 5M. until, for 
L, measurements of S/Al again agree with (5A4. Thus as \uh\ is increased from zero, the magnitude of the 
finite-size error SA4l — SM associated with a given choice of L, first increases from zero, reaches a maximum and 
then falls back to zero. 

As far as the number density itself is concerned, one finds |22| that on the line uh — 0, P converges rapidly to its 
limiting value with increasing L like pl — P ^ /^-(i-")/!^. Por nonzero \uh\, the dominant source of finite-size error is 
that described above. Thus for a given L, the finite-size error in the measured number density is minimised for state 
points on the line uh = 0, representing the analytical continuation of the coexistence curve io T > T^. Given that 
the critical isochore and the line uh = meet at the critical point and separate from one another only weakly as T 
increases from Tc (a feature that we have confirmed numerically), one can expect that the scaling field uh is generally 
very small on the critical isochore (see also Sec. 2) and hence that the system size dependence of pl{Pc, L) is less than 
it would be on any sub- or super-critical isochore within the critical region. 

To determine the finite-size critical isochore, the following procedure was employed. Taking a large, cubic, periodic 
simulation cell of linear size L — 17.5(t, runs were performed at the bulk critical point, for which p2| in reduced 
LJ units Tc = 1.1876(3), /ic/fcsTc — —2.778(2), pc = 0.3197(4). At this density, the system contains, on average, 
N — 1715 particles. Histogram reweighting |2^] was then employed to estimate for T = 1.0028Tc, the p value for 
which p = pc- A second simulation was then performed at this new state point, followed by a further extrapolation of 
the results to obtain pl{Pc, T) for the still higher temperature of T = l.OlSTc- This procedure was then iterated atotal 
of six times, giving isochoric data at temperatures T = T^, 1.0028rc, l.OlSTc, l.OSSTc, 1.087rc, 1.123Tc. Combining all 
six sets of simulation data via the multiple histogram reweighting scheme p6| , it was then possible to map the entire 
isochore in the range Tc < T < 1.123Tc. Fig. ^ shows the results, together with the estimate of the critical isochore 
for a smaller system of linear size L = lOfi, determined using an identical procedure. The differences between these 
two estimates of pl{Pc,T) (inset of fig. ^) are less than one part in 10*, confirming the expectation that finite-size 
effects are small. 
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FIG. 4. The measured isochore fiL{pc,T) for the two periodic system sizes L = 17.5a and L — 10a. The inset displays 
the difference between the two estimates, jj. is the absolute chemical potential (in units of ksT) subject to the convention of 
choosing the thermal wavelength A = 1 in the general definition (see Ref. jl^). Results are given in terms of LJ reduced units. 



C. Studies of the slit-pore geometry on the critical isochore 



Simulations of the LJ fluid conflned to a mesoscopic slit pore (cf. section IV A ) have been carried out for state 



points along the critical isochore of fig. IJ. Most of our studies are for a system of dimensions = Ly = 15a, L — 20a, 
although some results have also been obtained for a system of size — 65a, L — 20ct in order to gauge the magnitude 
of finite-size effects associated with the wa\\ area A = L^. We note that in the experiments of ref. |^ for SFg in CPG 
glass, the pore diameter is about 31 nm, corresponding to about 100 molecular diameters. 

Simulation runs comprised 10^ Monte Carlo sweeps for equilibration followed by 5 x 10^ sweeps for data collection. 
Measurements of the density profile were accumulated every 100 sweeps with each sweep involving (L^/rc)^ attempted 
particle transfers and (L^/rc)^ particle translations. Results have been obtained for both of the wall potentials given 
in equation (|4.2|), and for a range of values of the wall potential strength /. 



1. The critical point 

Density profiles p{z) at the bulk critical parameters fic, are shown in fig. ||(a) for the 4-10 potential and fig. ^(b) 
for the 3-9 potential. In both cases the bulk critical density pca^ — 0.3197 is denoted by a horizontal dotted line. 
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FIG. 5. (a) The measured density profiles p{z) (in units of a^) for the 4-10 potential at the estimated values of the bulk 
crit ical parameters fJ,c,Tc. The system size is = 15cr,L = 20a. Results for a number of wall strengths / are displayed (cf. 
eq (|4.2|)) (b) Same as (a) but for the 3-9 potential. In both cases the horizontal dotted line denotes the the critical density pc- 



It is instructive to compare and contrast the density profiles for both forms of wall potential. In all instances there 
is oscillatory structure close to the walls arising from excluded volume 'packing' effects. The number density at the 
wall (as measured e.g. by the height of the first peak) is principally controlled by the strength of the wall potential, i.e. 
by the value of /. Larger values of this parameter lead to a larger wall density and greater amplitude of oscillations 
in p{z). Further away from the walls, the packing effects gradually die out and p(z) varies smoothly with z. Within 
this smooth region, three main regimes of behaviour (common to both wall potentials) can be identified as / is varied. 
We address them in turn. 

For strongly attractive wall potentials (/a-g > 0.41,/4_io > 0.77), the density at the wall greatly exceeds the 
critical density and the profile is convex downwards with respect to the critical density. Thus the local density p{z) 
everywhere exceeds pc and decreases with increasing distance from the wall. At no point does it fall to pc- This latter 
feature is in stark contrast to the simulation results of Schoen et al |^J^ who investigated a very similar model to that 
described here. For certain fluid-wall interaction strengths, they reported a density profile that greatly exceeds pc near 
the walls, but dips below pc in the slit middle. Recently, however, this depletion feature has been demonstrated to be 
an artifact [ pTf , arising jointly from systematic errors in the simulation procedure used, and an incorrect designation 
of the critical point parameters |^ . 

For weakly attractive fluid-wall interactions (/a-g < 0.24,/4_io < 0.53) the density at the wall is less than the 
critical density and the profile is convex upwards with respect to the critical density. Thus the local density p{z) is 
everywhere less than pc and increases with the distance from the wall. At no point does it attain the critical density. 
Tests performed at subcritical temperatures indicate that for / values in this regime, the walls prefer the gas phase 
at coexistence. This is in contrast to the systems studied experimentally in ref. , where the liquid wets the walls 
for T <Tc and which, presumably therefore, correspond (in the language of the present model) to the large / regime. 

Turning now to intermediate wall strengths, the observed behaviour is somewhat subtle. Close to the wall, the 
packing-induced density oscillations span the critical density. Further away from the wall, the magnitude of the profile 
curvature is generally less than in the large or small / limits. Interestingly, there exist / values in this regime for 
which p{z) exceeds pc, but the profile is concave upwards with respect to the critical density. Thus p[z) exceeds pc 
and increases with increasing distance from the wall. We shall return to discuss this unexpected finding in section ^. 

Although both types of wall potential exhibit the same qualitative behaviour in the three regimes of / described 
above, differences are present in the detail. This is particularly true for small values of / as evidenced, for example, by 
a comparison of the profiles for /a-g = 0.18 and f^-io = 0.53. For these profiles the local densities in the slit middle 
are almost equal, but the wall density for the 4-10 potential is considerably greater th an tha t for the 3-9 potential. 
This difference reflects the relative range of the two potentials (as discussed in section IV A) and in particular, the 
fact that to obtain a given magnitude of wall potential at the slit middle (z = 10), one requires /4-10 = lO/3-g. 

We round off this subsection with some remarks concerning finite-size effects associated with the finite wall area 
A = L^. In a simulation, the infinite slit-pore limit L^j L 00 cannot be realised for all L values of interest because 
of bounds on the computationally accessible system sizes. Well away from criticality, this should elicit no grave 
concern because periodic boundary conditions in directions parallel to the walls provide a good approximation to the 
thermodynamic limit. A critical system, on the other hand is always 'aware' of its boundary conditions, by virtue 



14 



of its infinite correlation length. Changes in will therefore alter the effective range of correlations parallel to the 
walls, which might be expected to couple to those perpendicular to the walls in such as way as to affect the density 
profile p{z). To investigate this possibility, we have performed simulations in which we increased from the value 
Lx — 15(7 considered hitherto, to — 65(t. Owing to the high computational cost associated with such a large 
simulation cell, it was feasible to perform this comparison for only one type of wall potential, and we have chosen the 
4-10 form. The results are shown in fig. ^. Comparison with those of fig. ^(a) reveal that the increase in engenders 
only small changes in the form of the density profiles, the effect (such as there is) being greatest for intermediate wall 
strengths. On this basis it seems unlikely that the general scenario set out above would differ in the infinite slit limit. 
Nevertheless we feel that the role of wall area on confined critical systems is an issue that certainly merits a more 
systematic future investigation. 
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FIG. 6. As figure |5|(a), but for a system having Lx = G5a, L — 20a. 
2. Super-critical temperatures 

In this subsection we present results for the temperature dependence of the density profile on the critical isochore. 
Both the 3-9 and the 4-10 wall potentials have been studied in this regard. However, since it transpires that the 
qualitative form of the results are similar in both cases, we describe only those results pertaining to the 4-10 potential. 

Fig. shows the forms of p{z) at a selection of temperatures along the critical isochore of fig. ^, for the case 
fi-io = 0.89. This / value represents a strongly attractive wall potential, as evidenced by the high wall density. From 
the figure, one observes that at the critical point the local density p{z) is large compared to pc across the whole width 
of the slit. As the temperature is raised, however, the density in the slit middle decreases until, for T > 1.09rc, it 
reaches the bulk value pc- 
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FIG. 7. The measured density profiles p{z) (in units of (j^)for the 4-10 potential with /4-10 = 0.89 at a selection of 
temperatures along the critical isochore of figured. The critical density pc is denoted by a horizontal dotted line. 

The explanation of this behaviour is straightforward. In the vicinity of the critical point, the correlation length 
exceeds the slit width L and the density enhancement at the walls propagates across the whole slit, raising p(z) with 
respect to pc- In this regime one expects that for a sufficiently large slit width L, the density would decay to its 
critical value like z~^l^ ^ as is the case for critical adsorption at a single wall. Unfortunately our slit-pore is much too 
narrow for p(z) to reach the bulk value in the available range of z. Neither do we observe pure power law scaling for 
the variation that is present. This is because there exists no substantial region of z (away from the packing effects 
near the walls) which is not simultaneously influenced by the potentials of both walls. 

As the temperature is increased, there is a concomitant decrease in the bulk correlation length ^ t^^ which at 
some point becomes less than the slit width L. For distances from the wall that exceed ^, the decay of p(£) — p^. is 
expected to crossover to an exponential form, exp(— z/^). We indeed observe a rapid relaxation of pi^z^ to p^ for high 
temperatures (T > LOOTc), although for reasons similar to those described above, we have not been able to identify 
its character. 

In order to quantify the temperature dependence of p{z\ and to make contact with the experimental studies of 
refs and the t heore tical results of section |l|, we have obtained the temperature dependence of the adsorption, 
Tc, defined by Eq.(3.15) The form of Vc{j) is shown in fig. ^ for a representative selection of / values. The observed 
behaviour falls naturally into 3 regimes of / values, namely large, intermediate and small. 
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FIG. 8. The measured adsorption Tc = {p{z) — pc)dz ( in units of a^) for the 4-10 potential on the critical isochore, plotted 
as a function of the reduced temperature r = (T — Tc) /Tc. The slit width is L = 20a. Data is shown for a selection of / values, 
deriving from multi-histogram reweighting |26| of simulation data collected at six temperatures in the range Tc <T < 1.123Tc. 
For the two largest / values we also show fits to the data in the large r regime, of the form r(r) — a-f-br'^"'' with — u = —0.305. 
Statistical errors do not exceed the symbol sizes. 



For large / {fi-io > 0.77), fig. ^ shows that the adsorption increases monotonically as r is reduced to zero from 
above. Values of / in this range are believed to correspond to the situation studied experimentally in refs ||l|,||. As 
was the case for DFT results of Sec. [II, we find no evidence for the experimentally observed depletion phenomenon 



in which Tc{t) first rises to a peak as r decreases, and thereafter falls rapidly to negative values as the critical point 
is approached. 

We have attempted to analyse the form of Tc{t) for large / at temperature well above criticality. Assuming there 
exists a regime for which L ^ ^ 3> cr, (with a the particle diameter), one expects that Ti.{t) will exhibit the universal 
scaling behaviour of critical adsorption at a single wall, i.e. ^dT) ~ t^~^ . It is not clear, a-priori, that our rather 
narrow slit pore provides access to this regime. Nevertheless fig. || demonstrates that a fairly good fit to this form 
can be achieved in the high temperature regime. 

As / is reduced towards the intermediate regime, Tc{t) becomes progressively flatter until, for /4-10 ~ 0.65, the 
adsorption appears to exhibit little or no temperature dependence. This 'neutral wall' scenario presumably arises 
from a (near) cancellation of two competing factors associated with the wall potential. On the one hand there is 
the 'missing neighbours' effect (represented by the parameter c in equation (|3.4|)), whereby particles close to the wall 
have their potential energy raised relative to those in the bulk. On the other hand, there is the reduction in potential 
energy of particles arising from the attractive part of the wall potential (cf. the parameter in equation (3.4)). The 
cancellation of these two contributions effectively neutralises the influence of the wall on the fluid, and with it the 
temperature dependence of Tj,. This is also clearly visible in the corresponding density profiles (fig. M). 
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FIG. 9. The measured density profiles p{z) ( in units of cr ) for the 4-10 potential with /4-10 = 0.65 at a selection of 
temperatures along the critical isochore of figure 0. The bulk density pc is denoted by a horizontal dotted line. 



Turning finally to small / values, we find that the adsorption is negative at the critical point, reflecting the fact (cf. 
fig. I) that p{z) < Pc for all z. As r increases, however, Tc{t) is observed to increase monotonically with increasing 
T, although it at no point becomes positive. This increase in Tc is traceable to a progressive relaxation (away from 
the walls) of p{z) towards the bulk value pc. As in the case of large / discussed above, this effect has its origin in the 
decrease of the correlation length. 



D. A sub-critical isochore 

In the experiments of refs |l|,||, the adsorption was studied on sub-critical isochores having 0.995 < p/pc < 0.999. 
Unfortunately, in a simulation, it is not practicable to employ these same isochores because of the wholesale smearing- 
out of the critical region by finite-size effects. Specifically, in the FSS regime, one finds [H that L'^/''{p - Pc) ~ 
T{uhL'^^^'''), where T is some scaling function. In practical terms, this means that applying a given non-zero 
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bulk field to the fluid engenders a much smaller density change for a finite-size system than in the bulk limit. As a 
corollary, one finds that the typical scale of density fluctuations in a critical finite-sized system can be large. Thus for 
instance, for a cubic simulation cell of size L — 17.5a, the critical point density fluctuations extend from p « 0.15 to 
p « 0.5. Clearly therefore a fractional reduction in the density of 0.5% would have little discernible effect, either on 
bulk finite-size or adsorption phenomena. 

In view of this, we have chosen to study Tc{t) on the isochore having p = 0.22a^'^ = 0.68pc- Although this density 
is considerably smaller than pc, it nevertheless lies within the range of critical point density fluctuations in our finite- 
sized systems. One can therefore anticipate that if there is a critical depletion effect associated with negative values 
of the bulk field, then it should be visible on this isochore. 

The locus of the pa^ = 0.22 i socho re was determined in the temperature range Tc < T < 1.123Tc using the same 
procedures as outlined in section IV B. However, an additional complication in the present case was a greatly increased 
finite-size dependence of hl{p,TJ- Comparison of pl{p,T) for L — lOa and L = 17.5a revealed finite-size effects 
an orde r of m agnitude larger than those on the critical isochore (cf. inset of fig. ^). This feature was discussed in 
section IVB and is traceable to the differing magnitudes of the bulk field uh on each isochore. To ameliorate the 
problem, we determined p{p,T) using a very large cubic system of linear size L = AOa. Comparison of the results 
from L = 40(7 with those from a system of size L = 17.5a showed a discrepancy of less than two parts in 10^. 

The form of 1^c{t) for the 4-10 potential on the pa^ — 0.22 isochore is presented in fig. |l^. Data is shown for the 
same selection of / values given in fig. ||. Comparing these two figures, it is clear that the application of the negative 
bulk field reduces the adsorption to a considerably greater degree at large r than at small r. We attribute this to the 
large correlation length for small r, allowing the effect of the wall density (which is effectively pinned by the choice of 
/) to propagate across the slit, despite the action of the negative bulk field. The net effect is to increase the range of 
variation of Tc{t) for large / values but to reduce it for small / values. The same qualitative behaviour is observed 
for the 3-9 wall potential. However, in no instance do we see any sign of a rapid critical depletion of the pore density 
as T — > 0. 
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FIG. 10. The measured r dependence of the adsorption Fc = (pi^) ~ pc)dz ( in units of ct^) for the 4-10 wall potential on 
the pa^ = 0.22 isochore. Statistical errors do not exceed the symbol sizes. 



V. DISCUSSION 



The density functional (DFT) results and those of simulation have confirmed the heuristic arguments of Section 
II that there should be no critical depletion for a fluid in a single slit pore. In particular, for a fluid on its critical 
isochore, p — pc, and for strongly attractive walls the adsorption Tc increases monotonically as T is reduced towards 
Tc. There is no evidence of the rapid decrease of Fc near Tc which is the experimentally observed signature of critical 
depletion. Although the simulation results. Fig. 8, do not display the saturation of Tc that is clear in DFT results of 
Fig. 2 we believe this reflects the different values of L and the different temperature ranges used in the two treatments. 
In the DFT calculations, where L is very large, saturation does not set in until r ~ 10~^ and Tc is still decreasing 
quite rapidly with r for t > 10^"^. In the simulations, however, where L = 20a we would not expect saturation until 
T ^ To ~ 3 X 10"'^ and this is difficult to discern on the scale of Fig. 8. Such monotonic behaviour of Fc is consistent 
with the lattice gas results for bulk held = It differs from the Monte Carlo results of Schoen et.al. |4||^ but, 
as mentioned earlier, those simulations suffered from systematic errors p^ , p8| . The physical picture which emerges 
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from both our present treatments is one in which the density profiles show increasing enhancement, with respect to 
Pc in the central portion of the slit as the correlation length ^ increases and the decay at each wall becomes slower. 
When T ^ To, i.e. ^ ^ L, there is maximum enhancement and maximum adsorption. Increasing the strength of the 
wall-fluid potential (the parameter /) simply raises the overall level of the profiles and leads to an increase in Tc at a 
given value of r - see Fig. 8. 

For weakly attractive walls (small values of /) the simulation results yield negative values of the adsorption on the 
critical isochore, with Tc{t) monotonically increasing with increasing r. In this regime of / the wall-liquid interface 
would be dry for T < i.e. the walls prefer gas to liquid. At supercritical temperatures the density near the wall is 
lower than pc and the profile increases towards the centre of the slit. The closer one is to Tc the slower is the increase 
and the more negative is Fc- Similar behaviour of Fc is found for very weakly attractive walls (small Cu,) in our DFT 
calculations and was noticed earlier by Marini Bettolo Marconi By selecting an intermediate value of / it is 

possible to construct a nearly 'neutral wall' for which T(t) is small in magnitude and almost independent of r. In this 
case the density profiles are remarkably insensitive to r until very close to Tc where p{z) has the striking form shown in 
Fig. 9, i.e. it becomes concave upwards with respect to pc- More careful analysis reveals that for slightly greater values 
of / the profiles exhibit two symmetric maxima near the walls (but outside the regions of the oscillatory structure) 
for temperatures close to Tc. Similar, non-monotonic profiles were found recently for two-dimensional critical Ising 
films in the crossover regime between ordinary and normal transitions, i.e. for weak surface fields ||29| . Currently we 
are investigating whether the non-monotonic form of the profiles found in the present simulations of a Lennard- Jones 
fiuid arises from the same physical mechanism as in the Ising case. 

We were motivated to investigate sub-critical isochores because i) the sorption experiments on porous glasses were 
performed for r = {p — p^) / Pc = —0.001 and —0.005 and ii) the lattice gas studies of Ref. [^,0 suggested that depletion 
should occur on a path at fixed bulk field A/i = 2H < 0. Our present results for fluids show that no depletion occurs 
for the wide range of sub-critical isochores and pore-widths which we have investigated. The DFT results of Fig. 2 
imply that for the values of r pertaining to the experiments the behaviour of Fc(r) is changed very little from that 
on the critical isochore r = 0. Increasing | r \ does alter Fc(t) but this quantity remains monotonically increasing as 
T — > Tc. This is in keeping with the simulation results (with large /) of Fig. 10, although it should be noted that the 
simulations refer to r = —0.31 and to a much smaller value of L than in the DFT calculations or in the experiments. 
In summary, our explicit calculations appear to confirm the expectations of Section II that, for strongly attractive 
walls and parameter values of practical interest, the effective bulk field in the fluid is insufficient to drive Fc negative 
before ^ ~ L, i.e. we are in the regime <^ tq where depletion does not occur. Alternatively one can say that for the 
relevant values of r, the experimental pore widths L are not large enough for the density profile to relax from its high 
value near the walls to its bulk sub-critical value in the centre. Note that one cannot make |r| too large or one leaves 
the experimental critical region. The circumstances of the fluid are different from those of the lattice gas, where for 
'reasonable' choices of fixed bulk field H one has [|| tq '^th, the parameter equivalent to t^. 

We conclude that the adsorption calculated along critical and sub-critical isochores for a single slit pore does not 
exhibit the depletion observed in the experiments. What then can be the explanation of the experimental data? 
One might speculate that for some reason the density of the reference cell, which fixes the experimental isochore, 
was not constant but it is still difficult to see why this might mimic the fixed H scenario. What is more likely is 
that the observed depletion arises from the fact that the porous adsorbent is a complex solid material consisting of 
interconnected pores of various shapes and sizes whose morphology is poorly understood. Modelling such a material 
in terms of ideal (non-connected) slit pores is, of course, a gross over simplification. Recent theoretical and simulation 
studies pO[ show that the phase behaviour of fluids confined in such disordered media can be very different from that 
which occurs in a single pore but we are not aware of investigations of the adsorption along near-critical isochores. It 
should be possible to employ the model porous glasses produced in the quench simulations of Gelb and Gubbins , 
which mimic the spinodal decomposition process used to make real CFG glasses, in such an investigation. Were the 
adsorption to be vastly different from that found for the 'average' single pore this would be a striking demonstration 
of the importance of pore disorder and connectivity for critical phenomena. Note that further evidence that the single 
pore model is inappropriate to describe the experiments comes from the observation that critical depletion occurs for 
the colloidal graphite substrate Vulcan 3G for p/pc — 1-01 and 1.04 i.e. for super-critical isochores 0J^. Attempts to 
explain this observation within the single pore model require many further assumptions . 
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APPENDIX: A 



For the Landau model the integrals (|3.8| ) and (B.13) can be performed explicitly in the special case of A/i = (on 
the critical isochore pb = Pc)(see Ref. 11 5|): 



1/2 



(Al) 



where Vm = ip{L/2) — Pc)/pc is the reduced density at the midpoint of the slit, and D* = {DpV^)/kBTc. F{(j)\'nn) is 
the incomplete Jacobi elliptic integral of the first kind [p3] with arguments cos^ = rm/r^, where r^j = (p(0) — Pc)/Pc 
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is the reduced density at the wall, and m = (r^^^ + 2a* /b*) / {2r^^ + 2a* /b*). On the critical isochore pi, = pc the 
adsorption can also be evaluated ||l5| 



G = P\ 



2/3 



2D* 
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PcnL. 
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For A/i < 0, i.e. for p^ < pc, the relation (3^) between the order parameter at the midpoint r™ and the wall 
separation L can also be expressed in terms of the elliptic integral of the first kind. For T > Tc, a > one can also 
show 



L = 2p 
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2D* 
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{pq)-^^^F{cb\m), 



(A3) 



where the arguments of F{(f)\m) are given by tan(/)/2 — ^ 
real root of the cubic equation 



_ / q{rn,-rm) 
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and m = (l/4)M2):±(li:i^IlL. is the single 



+ + r{rl + 2a* /b*) + rf^ + {2a* /b*)r^ - (4/6*)A^* = 



(A4) 



and = {u- TjnY + "-^j 9^ = (^^ — fiY + where u = — (l/2)(ri + Tm) and = - - 2uri + 2a* /b*. 

Also for the case of Fisk-Widom free energy, the integrals in the formulae for L and Tq can be expressed in terms 
of elliptic integrals of the first kind for the special case of the critical isochore A/i = 0: 



L = 
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where y = r'^{{3a*/b*) + r^) and 



(A6) 



1 l.,,,2 1 1 



cos<^=(A^-- + -)/(A^ + ---) 
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and 



sin^a = 1/2- (3/4) 
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where 



A^ = (3(a76*) + 3rl)i/2, 



COS0 - (A^ - rl + r^)/(A2 + - rl) 



(AlO) 
(All) 



sin2a=l/2--(rf„/A2) 



(A12) 



Equivalent forms were derived in Ref. p5| | where they were used to study the behaviour of the adsorption in the 
crossover regime from the noncritical to the scaling region as a function of wall separation and surface fields. 
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